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ABSTRACT 



A single dimensional sinusoidal parameter perturbation adaptive 
control system is studied. The operation and signals at every point of 
the system are analyzed theoretically. Then, the system is simulated on 
analog computer starting with the simplest circuit and adding more com- 
ponents until a complete system is formed. The results of simulation 
study are compared with the theoretical analysis. Both agree in most 
cases. Finally, an analytical design of the system is attempted based 
on the simulation study and using a linearizing technique. It was simu- 
lated on analog computer and the performance was compared with the one 
predicted analytically. It does not follow exactly as predicted but the 
deviation is not so great. This method is useful in obtaining preliminary 
design which will be refined by experimental verifications. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Professor R. C. Dorf of the U. S. Naval Post- 
graduate School in this investigation. 
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1. Introduction 



The adaptive control system has received much attention recently. 

The principle of adaptive control Is to maintain a system always In an 
optimum operating condition whatever the affecting situation may be. There 
are two kinds of adaptive control systems; the signal adaptive system and 
the system adaptive system. The signal adaptive system is a system which 
can analyze incoming signals and automatically adjust the system to give 
an optimum output. The system adaptive system is a system which can mea- 
sure its own characteristics and then change its structure or design so as 
to continuously maintain a desired performance. This paper deals with only 
the latter. The system adaptive system should have the following two es- 
sential features to achieve its goal; 

(1) The measurement of its own characteristics. 

(2) The generation of proper corrective signal for the control- 
led system according to the measurement. In the measurement of the 
characteristics of the system, the optimum operating condition or figure 
of merit is defined. The figure of merit or index of performance (IP.) is 
an indication of the characteristic or quality of a system. The IP may be 
a physically meaningful value such as heat produced per gallon of fuel or 
theoretical quantity such as RMS error or mean squared error. 

The sinusoidal parameter perturbation adaptive control system is a 
method of system adaptive control and in this paper the operation and 
signals at every point of this system, are analyzed theoretically. Then, 
the system is simulated on analog computer starting from the simplest cir- 
cuit. Adding more components one by one, the effect of each component is 
observed until the complete system is formed. The results of these simula- 
tion study are compared with the theoretical analysis. Both agree in most 
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cases. Finally, a design of the system Is attempted analytically based on 
the simulation study and using a linearizing technique. It was also simu- 
lated on the analog and the digital computer and the actual performance of 
the system Is compared with the predicted performance from the analytical 
design. The results show that the actual performance does not exactly agree 
with the performance predicted analytically but It follows reasonably close 
pattern predicted analytically. This Is expected because the analytical 
method Involves approximation. This method Is useful to obtain a prelimin- 
ary design which will be refined by experimental verification and gives a 
direction of Improving the first trial design. 
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2. The Sinusoidal Parameter Perturbation Adaptive Control System. 

A single dimensional sinusoidal perturbation system is as Fig. 1. 

The model is an optimum operating model. When the plant is operating 
exactly as the model there will be no error output. But when a parameter 
in the plant is changed to some other from the optimum parameter value due 
to environment change, there will be an error output. By squaring the er- 
ror we get a non-negative function of error. As an index of performance 
we may adopt average of error squared as: 

i r T 2 

IP = i \ e dt (2-1) 

J 0 

This quantity is a function of the characteristics of input to the plant 
and also of the characteristic of plant parameter drift caused by environ- 
ment change. Fig. 2 shows the contours of IP(X,Y) for two dimensional 
case where IP is a function of two parameters X and Y. The minimum IP 
occurs at X=X Q and Y=Y If we fix the value of Y at Y^ and vary only X, 
the IP variation with X becomes as Fig. 3, the single dimensional case. 
Generally the IP with one parameter variable tends to be a parabolic func- 
tion of the parameter and can be approximated near the optimum parameter 
value as; 

IP(X) = ai (X-X 0 ) 2 + a 2 (2-2) 

where X^ is optimum parameter value and a^ and a^ 
are arbitrary constants which change the shape 
of IP curve such that it represents closely the 
actual IP variation of any plant. 

In some plants a parameter has a slight or no change with respect to the 
environment change and such parameter is fixed at the optimum value and the 
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Fig. 1. Single Dimensional Sinusoidal Parameter Perturbation Adaptive 
Control System. 
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number of parameters to be controlled is reduced. 

Referring to Fig. 3, suppose we perturb the parameter X sinusoidally 
by a perturbation signal, C p = sin oO^t. If the actual plant para- 
meter X is at the optimum parameter value, X^, there will be no IP varia- 
tion output at frequency U) y If the X value is lower than X Q , there will 
be an IP variation output at frequency CO^ and the output signal will be 
180® out of phase with respect to the input perturbation signal as seen in 
the Fig., on the contrary, if the X is higher than the X^, the output is in 
phase with the input perturbation signal. It is also noted that the ampli- 
tude of the output is proprotional to |x-X^|, therefore the output is consid- 
ered to be a carrier, modulated by the variation of |x-Xq| through the IP 
characteristic. The frequency of carrier is CO ^ and it carries the informa- 
tion about the polarity of the X-X^ by its phase difference with the input 
perturbation signal. When this output is demodulated by multipling with 
the original perturbation signal, it will produce a D.C. signal with proper 
sense to correct the plant parameter deviation. 

The bandpass filter eliminates unwanted signal components contained in 
the modulated signal and also affects the response of the adaptive loop to 
a step change of signal. 

The detector is a multiplier which multiplies the modulated signal with 
original perturbation signal and produces positive or negative correction 
signal. 

The filter after the detector is a smoothing filter which filters out 
high frequency components in the correction signal and acts as compensator 
to stabilize the adaptive loop. 

The integrator is a D.C. motor in practical operation and it integrates 
correction signal and drives the plant parameter control to the optimum value. 
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If we want to make more than one parameter adaptive, separate loops 
for each different parameter must be established and use different perturba- 
tion frequency for each parameter. The same error measure can be used for 
each loop if it is suitable, but for each parameter, other elements must be 
duplicated. 

Sometimes, a plant has several local minimum in its IP characteristic. 
In that case, a special measure is taken to find the true minimum. But 
generally, all plants have a single minimum around the optimum parameter 
value. The sinusoidal perturbation method can also be applicable to an 
adaptive system which searches for a maximum IP value. 
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3. The Analysis of Signals in the Adaptive System 
3.1 Block diagram of an adaptive system 

To set up a block diagram of a single dimensional parameter perturba- 
tion adaptive control system, the circuit models of each component are 
established. Referring to Fig. 1, the parameter perturbation signal per- 
turbs a plant parameter and the parameter itself also is changing its value 
due to environment change and parameter correction signal coming through 
adaptive loop. This parameter perturbation signal and parameter drift 
cause the error output while the plant and model are operating with input, 
R. It is necessary to find out an equivalent transfer function block be- 
tween the parameter perturbation input and the error output. The para- 
meter drift is also an input for the error output and the input is same as 
perturbation input. We may designate this equivalent transfer function as 
error measurement transfer function and express as 



K 



m 



G (s) = - 

m 1 + £» s 



(3-1) 



K and are decided by the signal transfer characteristic between 

perturbation input and error output. To find out the transfer character- 
istic, phase lag and attenuation introduced to this error output with res- 
pect to the original perturbation and parameter drift are observed. Al- 
though the perturbation and drift do not come into the plant through 
ordinary input point but through parameter itself, these signals encounter 
the same phase shift as a signal of the same frequency coming into the 
plant through the input of the plant, encounteres. This will be confirmed 
by actual simulation study in Section 4.1. The attenuation factor of this 
error signal is not the same as the attenuation factor of the signal of same 
frequency coming through plant input. The attenuation factor to the error 
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signal depends on the frequency and magnitude of the plant input, R and 
the characteristic of the plant. It is also governed by its original 
signal frequency. If we fix the amplitude and frequency of the input to 
the plant, the attenuation factor is a function of its original perturba- 
tion or drift signal frequency assuming the plant characteristic are con- 
stant. 

From the above discussion, the^ is chosen such that the equivalent 

tn 

transfer function gives the same phase lag as the plant introduces to the 

same signal coming through the plant input. 

Next, the K is decided by actual measurement of attenuation, 
m 

In measuring the attenuation, the input to the plant is kept constant and 
perturb or drift the parameter by some signal and measure the error out- 
put, which gives the attenuation factor between parameter drift input and 
the error output. Let the measured attenuation be A . Then 



A 

m 



K 



m 



a 



.r 2 

m 



) 



m 



(3-2) 



where CO is a parameter drift frequency used 
m 

for measuring A . 

m 

From equation (3-2) K can be calculated. 

m 

2 

The IP curve is represented by IP(X) = a^CX-X^) + a^ as explained in 
Section 2, A multiplier is used as squarer and a^ is the coefficient of 
the squarer. After the squarer, is added. A proper choice of a^ and 
gives an approximate representation of any single dimensional IP curve. 

The bandpass filter has the following form, 

co 1 s 

G.(s) = 2 2 (3-3) 

1 S + 2/«^S + U) 

The detector M2 is a multiplier. 
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The smoothing filter G 2 (s) takes various forms depending on the re- 
quirement of compensation of the adaptive loop but generally it is a phase 
lag filter. 

The reference signal filter is an arbitrary filter chosen such that 
phase shift at frequency faj ^ due to this filter is equal to the phase 
shift introduced to the perturbation signal between the perturbation in- 
put and detector through the adaptive loop. 

The integrator is usually a D. C. motor and its transfer function is 
designated as G^(s). 

The complete block diagram of a single dimensional parameter perturba- 
tion adaptive control system becomes as Fig. 4. 
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Fig. 4 Block diagram of a single dimensional sinusoidal 
parameter perturbation adaptive control system. 



3.2 Modulation 



\ 



The input signals through the parameter are the perturbation signal 
and parameter drift signal due to the variation of the parameter value. 
The instant error deviation of parameter, Y^(t) can be expressed as 

Yj(t) = X(t ) + C p - X Q (3-4) 

where X(t) = instantaneous parameter drift 
Xq = optimum parameter value 
C p = parameter perturbation signal 

The analysis of signals will be made for two cases of parameter change; 
sinusoidal drift of parameter and step drift of parameter. 

(a) Sinusoidal drift of parameter 

The parameter drifts sinusoidally around the optimum parameter value 
with very low frequency. C = C,sin£c\t is used as perturbation signal 

p 1 1 

throughout the analysis. The instantaneous parameter drift is 

X(t) = X Q + X x sina> m t (3-5) 

CO = parameter drift frequency and is much lower 
m 

frequency than perturbation frequency, u)^ 

Then Y^(t) becomes as below using equation (3-4), 

Y, (t) = C-sin t + X, sin^y t (3-6) 

1 1 1 1 m 

Referring to the Fig. 4, Y^(t) goes through G m (s) and comes out as 
Y 2 (t), 

Y~(t) = e(t) = A C.sinC^y.t - 0 ) + A X..sin(£tl t - 9 ) (3-7) 

Z pi 1 p ml mm 

where A , A are associated attenuation factors, and 

p* m 

9 , 0 are associated phase shifts. 

p m 
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Signals at each point in the loop are labeled in the Fig. 4. 
Following the loop, 



V'> 



- a l [ V 2< e >] 



Y 4 (t) = a 2 + Y 3 (t) 



Calculating above equations, Y^Ct) becomes 



V fc) 



a 2 + a i ( ' 



, 2„2 
A C . 
P 1 



. 2,2 ,2 2 

A X, a, A C, 0 

m 1 . 1 p 1 cos2 

- + - 2 ' ' 2 
2 2 

a i m 1 cos2(&) t - 9 ) 

2 m m' 



(3-8) 

(3-9) 



<“Y - V 

(3-10) 



+ 2a,A A C,X,sin(a) t - 9 )sin(^>,t - 9 ) 

• 1 m p 1 1 m m 1 p 

The last term of Y,(t), 2a, A A C,X,sin(£t> t - 9 )sin(«/.t - 9 ) is the 

4 lmpll m m 1 p 

modulated signal term where the perturbation signal is amplitude modulated 
by parameter drift signal. This term carries not only the information of 
magnitude of drift signal but also the information about the actual operat- 
ing point of parameter, in other words, the actual parameter is higher or 
lower than the optimum parameter value by the following form in an actual 
signal. 

When X(t) - X A ) 0 or X-sin(&; t - 9 ) > 0, this term is 
0 1 mm 

2a A A C I X. sin( Ut t -9)1 sin(^t - 9 ) (3-11) 

1 m p l » l m mi 1 p 

When X(t) - X A 0 or X-Sin(&> t - 9 )/ 0, this term is 

2a A A C I X- sin ( Co t - 0 ) I sin(&;,t - 9 - n) (3-12) 

1 m p 1 | 1 v m nr I v 1 p 7 

This is explained in Section 2 and shown in Fig. 3. 
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(b) Step drift of parameter 

The parameter value drifts in step from optimum parameter value to 
some value X, then by equation (3-4) , the instantaneous deviation of para- 
meter becomes 

Y.(t) * X + C - * (X - X J + C.sin&At (3-13) 

1 p U U 1 1 

Going through the equivalent transfer block, G (s) and comes out as 

m 

Y 2 (t) = e(t) = A d (X - X Q ) + A CjSin^t - 0 ) (3-14) 



where A,, A are associated attenuation factors and 
d* P 

0 is the phase shift. 

P 

Calculating the same way as in case (a) using equations (3-8) and (3-9), 
Y^(t) comes out as 



a 2 c 2 

V*> * a 2 + a i ( A d< x - V 2 + ~h 



. 2 2 
a, A C, 






co scO t 
1 



+ 2a 1 A d A p C 1 (X - X Q )sin(&; i t - 0 ) (3-15) 

In Y^(t) the last term, 2a^A d A p C^(X - Xg)sin(tt/jt - 0 p )is an amplitude modulated 

signal term where the carrier is the perturbation signal, sin {CO, t - 0 ). 

1 P 

This term carries the information about the amount of parameter drift (X-X ) 
and furthermore the direction of parameter drift from the optimum parameter 
value. This term takes following form for actual signals: 

When X - X^ > 0 or the actual parameter is higher than the optimum 
parameter value, this term becomes 

2a i A d A P c i i x - x 0 1 sin <“’i t • y < 3 - i6 > 

When X - X^ 0 or the actual parameter is lower than X^, this becomes 
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X 



(3-17) 



2a A A C. 
1 d p 1 



as explained in Section 



- X Q | s in ( CO^ t - e ? - it) 
2 and shown in Fig. 3. 
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3,3 Filtering 



The bandpass filter G^(s) attenuates considerably all unnecessary 
signals in Y^(t) and only passes the components that contain information 
about the parameter drift* without much attenuation. 

(a) Taking the case of sinusoidal parameter drift* the Y^(t)* equation 
(3-10) in Section 3.2(a), is rewritten as 



Y 4 (t > = D 0 ' D 1 cos2 (^ 1 t - G p ) = D 2 cos2(^ m t - 9 m ) 



+ D 0 sin(c«; t - 0 )sin(t<J,t - 6 ) 
3mm Ip 



(3-18) 



where 



a. 



= 



1 2 2 2 2 

a n + -^-(A C + A XT) 
pi ml 



0 2 2 



D, = 



A 2 2 
a A C. 
1 P 1 



D„ = 



2 2 
a A X, 
lml 



D 0 = 2a- A A C,X- 
3 1 m p 1 1 

Expanding the last term of equation (3-18) 5 



Y 4 (t) = D Q - D 1 cos2(io 1 t - 9 p ) - D 2 cos2( a ; m t - 0J 

s D 



3 cos f 

2 L 

) 



(CO. - (0 ) t - 9 + 0 1 - COS \(CU. + UJ )t - 0 

% 1 nr p ml 2 pirn 7 p 



(3-19) 



The Y^(t) contains the following frequency signal components 5 

D.C. , 2 U> ,2 CO , CO - 0) „ CO + lO . 

* m* 1 1 m 1 m 
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The necessary intelligence is contained in CO. - UJ and CO + U) com- 

I m I tn 

ponents. Therefore the bandpass filter is designed to pass these two 
components without attenuation and suppress other components. Thus the 
center frequency of the bandpass filter should be CO ^ and this filter does 
not introduce any phase shift to the signal of Assuming that the 

phase angle characteristic is linear near the center frequency, which is 
generally true for an actual filter, the signal of ( CO ^ - U gets + A 9 
phase shift and the signal of ( CO ^ U> m ) gets - A 0 phase shift with 
respect to the signal of Then Y^(t) passing the filter, becomes as 

• Vo - A lV° s2 ( “V - 0 p - V - V 2 - e „ • V 

+ _V3_co S (" ,to )t - 9 p + 0 m+ A e] 

- —2 r - co ! (<"i + «V* - 8 p - 9 „ - * 0 ) < 3 ' 2 °) 

where A„, A,, A.,, A., and A. are attenuation factors associated eith each 
0 1 l 3 4 

signal and 

0j, 02 and A 0 are additional phase shifts introduced to each signal. 
Usually in a band pass filter A^ and A^ are the same attenuation as the 
center frequency signal suffers because as assumed in Section 3.2. 

A^, A^ and A 2 should be much smaller than A^ for our purpose. Combining 
the last two terms of equation (3-20), Y^(t) becomes 

V 5 (t) * Vo • Vl cos2 <‘ u l t - 9 p - V • A 2 V os2 <"n, t - \ ■ V 

-f A,D.sin( CO , t - 0 ) sin ( to t - 9 - A 9) (3-21) 

3 3 1 p mm ' 
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If the modulated signal term in Y^(t), equation (3-18) is compared with 
the modulated signal in Y^(t), equation (3-21), it is noted that the band- 
pass filter introduced a phase lag of A 9 only to the modulating signal 
component, (sin - 9^) as far as the modulating signal term is concerned. 

(b) In step parameter drift case, Y^(t), equation (3-15) is rewritten as 



Y 4 (t) = D $ - D 6 cos2 Ct ) 1 1 + D y (X - X^sin^t - 0 ) 



( 3 - 22 ) 



where + a^ 





D, = 2a A A C. 
7 1 d p 1 



Passing the bandpass filter, Y^(t) becomes 



Y $ (t) = A q D 5 - A 1 D & cos(2 UJ 1 t - 0 2 ) + A 3 D ? (X - X 0 )sin(<^ 1 t - 9 p ) (3-23) 
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3.4 Demodulation 



The modulated signal Y<.(t) reaches the detector where Y^(t) is 
multiplied by the detection reference signal and demodulation is accomplish- 
ed. The detection reference signal should be in the same phase as the car- 
rier signal to get maximum output at the demodulation. For this purpose, 
detection reference signal phase lag filter G^(s) is used to give the refer- 
ence signal an exactly same phase lag as the phase lag introduced to the 
carrier signal between the perturbation input at the plant and the detector; 
in our case which is 0^. Thus the detection reference signal has the follow- 
ing form, 

Y 6 (t) = BC lS in( (V 1 1 - 0 ) (3-24) 

where B = attenuation factor of G^(s) 

0 = phase lag introduced by G_(t) 

P 3 

(a) For the sinusoidal parameter drift case, Y,.(t) equation (3-21) is 
multiplied by Y^(t) equation (3-24) and we get 

Y y (t) = A Q D 0 BC 1 sin( IV ^ - 0 ) - B C 1 cos2(^ 1 t - 0 p - 0 p 

sin( Co. t - 0 ) - A-D.B C,cos2(o> t - 0 - 0_)sin(a4t - 0 ) 

1 p mmz 1 p 

+ A_D_B C.sin 2 (£c> 1 t - 0 )sin(a; t - 0 • 40) (3-25) 

331 lp mm ' 

In Y y (t) only a D.C. component is integrated by the integrator following and 
develops a parameter correction signal. Only the last term of Y y (t) contains 
a D. C. component which is varying slowing in accordance with the modulat- 
ing signal, which is the parameter drift. All other terms do not contribute 
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to the parameter correction because they are comparatively high frequency 
sinusoids and furthermore, A^, and A^ are very small due to the band- 
pass characteristic of G^(s) and they will be further attenuated by the 
smoothing filter following. If only the last term of Y_(t) is expanded 
it becomes 

A T) R C ^1 

3 3 lsin ( to t - G - 4 6) - -sin( to t - 9 -AS) 

2 mm 2 v m m 

cos2 (CO - 9^) (3-26) 

The first term of (3-26) is considered as a slowly varying D. C. which 
corrects the parameter drift. Thus, 

A D B C 

Y (t) = sin( COt - 9 - A 9) (3-27) 

/ d ♦ c i m m 

Putting back the original value for and taking into account the phase 
shift between Y^(t) and Y^(t), 

Y Q (t) = a. A A A.B C^X,sin( t - 0 - A 6 - 0.) (3=28) 

y 1 m p 4 11 m m J 

where 0^ is the phase shift in G^Cs) anc * 

A. is the attenuation factor associated 

4 

The parameter correction signal Y^(t) has the same form as the parameter 
drift and has some additional phase lag compared to the parameter drift. 
This correction signal will be fed back negatively to parameter adjustment 
and corrects the drifted parameter to the optimum parameter value, 

(b) For the step parameter drift case s Y^(t) equation (3*23) is multiplied 
by Y^(t) equation (3-24) and we get 

Y_ (t ) = A D C B C, sin( tV , t - 0 ) - A-D.B C-eos^ft* t - 0 O ) 

7 Ool 1 p 101 I A 
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(3-29) 



+ A D B C.(X - X.) sin 2 ( u>,t - 9 ) 

J / 1 u 1 p 

In Y^(t) only the last term develops a I). C. component, which will be 
integrated and drive the parameter adjustment. Other terms do not contri- 
bute to parameter correction because they are comparatively high frequency 
sinusoids and and are very small due to the bandpass filter attenua- 
tion and will be further attenuated by the smoothing filter. If only the 
last term of Y^(t) equation (3-29) is expanded it becomes as, 

— A 3 D y B C^X-Xq) - 4“ A 3 D 7 B C 1 (X-X 0 )cos2(^ 1 t - 0 ) (3-30) 

Putting back the original value for D^, 

Y 7d.c. (t) * a l A d A p A 3 BC ? (X - V < 3 ' 3l > 

Thus the rate at which the parameter drift is corrected, is proportional to 
the parameter error (X-X^), the perturbation signal amplitude C^, the con- 
stant a^ associated with the IP curve (the sensitivity of the IP variation 
to changes in parameter value) and adaptive loop gain. This open loop 
analysis is not exactly correct when the loop is closed, but it shows a 
possible signal development through the adaptive loop. 

(c) Next, the effect of phase difference between the carrier and the detec- 
tion reference signal Y.(t), in the demodulation, is considered. To see 

o 

more clearly the process of demodulation and the effect of the phase dif- 
ference between carrier and the detection reference signal, graphical multi- 
plication of actual signals involved, is performed. Taking the case of 
step parameter drift, the actual modulated signal in Y^(t) has the following 
forms as explained in 3.2(b) and expressed in equations (3-16) and (3-17), 
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when X - 0 



(3-32) 



A X - X I sin(£c/t-0 -n) when X - 0 

I u | 1 p 0 ^ 



(3-33) 



where A is a constant. 



The carrier in Y,-(t) is C^ 5 (t) = sin( CO^t - 0^) when X-X Q > 0 and -C^ 5 (t) = 
sin(^-t - 0 - it) when X-X^ < 0. 

1 p ' 0 ^ 

The detection reference signal is Y^(t) = B sin( CO t - 0^). 

If 0^ ^ 0^ Y^(t) l ea< *s or lags the carrier C^(t) by some angle. The 



graphical multiplication is performed for the following three cases* 



The coefficients, of each term is a constant* therefore only a sine wave of 
unit maximum amplitude represents each wave. 

Fig. 5 shows the demodulated output by dotted line for the above three 
cases. It shows clearly the phase sensitive nature of perturbation signal; 



correct direction of parameter correction in accordance with the actual 
parameter error* X-X^^ 0 or X-X^ < 0. 

In case (1), we get maximum net output of correction signal when inte- 
grated as seen in Fig. 5(1). For the case (2)* we get less net output of 
correction signal when integrated. This eventually ends up with less adap- 
tive loop gain. The case (3) has the same output as in the case (2). 

If we consider the time lag of starting to build up correction signal 



(1) Y A (t) is in phase with the carrier. 



(2) Yg(t) leads the carrier by 30°. 

(3) Yg(t) lags the carrier by 30°. 
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output at the integrator, the case (3) has the largest time lag of the 

three cases, and the case (1) is the second and the case (2) is the least. 

But this time lag is negligible because the carrier frequency is very high 

compared to the modulating frequency. 

From this analysis, we can see that when the phase difference between 

Y c and Y, is exactly 90°, the signal output will be the lowest and when the 
o b 

phase difference is more than 90° the demodulated signal polarity will be 
reversed and the magnitude will increase as the phase difference approaches 
180°. Therefore, when the carrier phase lag is expected to be closer to 180°, 
it is better to remove the reference signal phase lag filter and change the 
polarity of the demodulated output to have correct feedback. 
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4. The Simulation Study of a Sinusoidal Perturbation System 

4.1 IP. variation investigation 

To investigate the actual variation of the IP. with respect to the 
parameter drift and perturbation of the plant 5 the following model and 
plant were simulated by an analog computer with the adaptive loop opened. 
The analog computer used throught the Investigation is the Bonner 3100. 



Model; ^ ^(s) 



K 

S(S + x 0 ) 



Plant ; G p (a) - s(s + x) 



The pole of the plant is a variable parameter x 9 and the pole of the model 
is fixed at 3, which is the optimum parameter value in this case and gives 
"jf of 0.5. The pole of the plant can vary in the range that gives jf of 
plant from 0.1 to 1. with K » 9. The block diagram of the simulation is as 
Fig. 6 and the analog computer circuit is shown in Fig. 7. 

(a) The input to the plant and model is a square wave of 0«1 cps, 

A plant parameter, the pole of the plant is drifted by a triangular wave 
of 0.01 cps. while the model and plant are operating. The output of model 
and plant are as Fig. 8(a) and (b) respectively. The output wave shape 
shows the change of jf of the plant due to the parameter drift. 

(b) The error output is recorded as Fig. 8(c) and the pole drift 
is as Fig. 8(d). 

Next, the input to plant and model is sine wave of 0.3 cps. and the 
error output is also recorded as Fig. 9(a). The Fig. 9(b) is the pole drift 
as before and (c) is the error squared output. 

The above results shows that error is produced only when the output 
of plant is in transient and the plant response Is not the same as the 
model response due to the parameter drift. Also it shows that the magnitude 
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Figo 6 The model and plant set-up for the investigation of XP 0 
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Fig. 7 The analog computer circuit of the model c nd plant set-up for the 
investigation of IF. 




(c) 



(d) 




Pig. 8. The output of model, plant and error.with plant pole drifted. 

(a) model, 5 

(b) plant 

( c ) error 

(d) plant pole drift 



l 
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Pig. 9. Error and error squared output with sine wave input to plant. 

i a) error 

b) plant pole drift 
(c) error squared 



l 
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of error is proportional to the parameter error, the magnitude of para- 
meter drift. 

The magnitude of error is bigger when the parameter drift is in the 
direction of decreasing plant with the same amount of parameter drift. 
Therefore, it is expected that the magnitude of parameter correction signal 
will be greater for the parameter drift to a decreasing J* than a increasing 
j’ , with the same amount of: drift. 

(c) Error output is observed with only perturbing the pole of 
plant by a sinusoidal perturbation signal. The perturbation frequency is 
much higher than a drift frequency; here a sine wave of 5 cps is used. 

(1) The magnitude of error output due to perturbation signal 
only. The perturbation wave appears in the error output, but the magnitude 
is much reduced due to the severe attenuation. Five cps is a high frequency 
for this plant. Thus the parameter perturbation signal does not appear con- 
spicuously in the plant output even with a big amplitude perturbation signal 
as far as the frequency of perturbation is high. If we want to eliminate 
the effect of perturbation on plant output completely, the parameter of 
model can be perturbed instead of plant parameter. Fig, 10(a) is the per- 
turbation wave appearing on the output of the plant. The perturbation 
waves appear conspicuously on the output of the plant only when the ratio of 
the amplitude of the plant input to the amplitude of the perturbation signal 
is small. 

(2) Phase shift of perturbation signal through plant. 

The perturbation wave appearing on the output of plant lags the original 
perturbation signal by some angle. This phase lag is introduced by the 
plant and depends on the characteristic of plant and signal frequency. 

Fig. 10(c) and (d) are the original perturbation signal and the perturbation 
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(a) Perturbation wave 
appearing on plant 
output* 



(b) Original perturbation 
signal. 






(c) Original perturbation. (e) Input to plant 

(d) perturbation wave on output (f) Output of plant 



Pig. 10. Perturbation waves on the output of plant. 
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wave appearing on the plant output respectively. The same marks on the 
waves are the corresponding points on both waves. These corresponding 
points are found by starting and stopping the signal abruptly. According 
to the measurement of phase difference between the two waves, perturbation 
wave appearing on the output lags the original signal by 165°. To compare 
this phase lag with the phase lag introduced by plant to a signal coming 
into the plant through the plant input, the same signal used as perturbation 
is fed into the plant. Fig. 10(e) and (f) are the input to plant and the 
output of plant. The same marks on the signals are the corresponding points 
on both waves which is found as before. According to the measurement, the 
output waves lags input by 165° approximately, which is the same phase lag 
introduced by plant to the parameter perturbation signal. From this re- 
sult it is concluded that the phase lag introduced by the plant to the 
parameter perturbation signal fed through the pole of plant is the same 
phase lag introduced by plant to a signal of same frequency fed through the 
plant input. 
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4.2 The Simplest Adaptive System Simulation (Adaptive circuit 1) 

The complete adaptive control system is as Fig. 1 and the block dia- 
gram is as Fig. 4. To investigate the function of each component and the 
effect of change in components on the over-all system operation, the simula- 
tion is started with the simplest combination of components for a adaptive 
system. Then more components are added to improve the system operation 
until the complete adaptive system is formed. We designate successive cir- 
cuits by number from the simplest to the complete circuit in simulation 
study. The first circuit is designated as adaptive circuit 1, which is the 
simplest and is shown in Fig. 11(a). The analog computer set-up is as Fig. 
11 (b). It simulates starting from the error output point in the complete 
system. In Section 3.2, signal analysis, the instant deviation of parameter 
Yj(t) is expressed as equation (3-4), 

(t) = X(t) + C p - X Q (4-1) 

For the analysis the optimum parameter value X^ is arbitrary and can be 
set as zero. Y^(t) causes an error output and the error output Y 2 (t) has 
severe attenuation and phase shift and comes out as equation (3-6) for a 
sinusoidal parameter drift. But in the simulation study, ignore the at- 
tenuation and phase shift and assume the error output Y 2 (t) as, 

Y 2 (t) = X(t) + C p (4-2) 

This error output is going to be the input in the Adaptive circuit 1 because 
the circuit simulates from error output point. In an actual system, the 
parameter correction signal corrects plant parameter and reduces the error 
output to zero but in this simulation the correction signal directly reduces 
the Y 2 (t), which is equivalent to correcting the plant parameter. 
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(a) block diagram 




Pig. 11. Adaptive circuit 1. 
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In the circuit, and M 2 are a squarer and a detector, both of which 

are electronic multipliers. The zero D. C. level of electronic multiplier 

is not constant, therefore a pure integrator can not be put in right after 

the multiplier because the integrator integrates the drifting D. C. and 

K 4 

developes erratic correction signal. A pseudo integrator - — — is used 

v S + 0 . 1 

4 

for the simulation of the pure integrator — . The parameter perturba- 

b 

tion and parameter drift signals used are as follows: 



C = C.sin^.t 

p 1 1 


= 1 v. 


*>1 


= 5 cps. 


X(t) = X.sinto t 
1 m 


Xj = 8 v. 


to 

m 


=0.1 cps. 



The adaptive loop is open at Yg feedback point and open loop observation is 
made . 

(a) The wave shapes of signals at every point of loop are shown 
in Fig. 12. The (a) in the Fig. is the error signal in which perturbation 
signal is superimposed on the parameter drift signal. (b) is the output of 
squarer where the modulation is accomplished. (c) is the output of detector; 
the demodulated output has positive or negative D.C. component in accordance 
with the parameter drift direction. This is integrated and comes out as (d), 
which is the parameter correction signal. If the parameter correction sig- 
nal is compared with the parameter drift signal in (a), it has the same shape 
as the parameter drift, therefore, when the correction signal is negatively 
fed back it will correct the parameter drift. But the correction signal lags 
the drift signal by 60° approximately as seen in (a) and (d). The marked 
points are the corresponding point on both waves. This phase lag is caused 
by the integrator as seen between (c) and (d). There is no phase lag be- 
tween Yj and Y^ as seen in (a), (b) and (c). The perturbation signal ap- 
pears in the correction signal Yg, which will be eliminated by adding a 
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(a) Error, Y 2 (t) 
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(b) Error squared 
Y 3 (t) 




time 



(c) Demodulated outp ut Q 



Y ? (t) 
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(d) Correction signal 0 



T 9 (t) 




Pig. 12. Wave shapes of signals in Adaptive circuit 1. 
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smoothing filter, 

A square wave also can be used as a parameter perturbation signal. 

The same results as obtained by sinusoidal perturbation, are obtained 
using square wave perturbation. The square wave perturbation wave ap- 
pearing in the correction signal looks like a sinusoidal wave. Fig. 12-1 
shows the results, 

(b) When the frequency of perturbation is increased, it is at- 

tenuated more and it does not appear in the correction signal; in this 
case there is no perturbation signal in the correction signal when 15 cps. 
is used for ^he am Pl* tu< * e °f the correction signal is not affected 

by the change of perturbation signal frequency. This means the perturba- 
tion frequency does not affect the loop gain in this simple circuit, how- 
ever, in an actual circuit the frequency of perturbation will affect the 
amount of attenuation due to circuit components and affect the loop gain. 

Fig. 13 shows the correction signal with different perturbation frequencies. 

(c) Next, the perturbation frequency is fixed at 5 cps. and the 
frequency of parameter drift is varied. When the drift frequency increases, 
the amplitude of correction signal is decreased. Fig. 14 is the amplitude 
of correction signal vs. drift frequency graph. The decrease of amplitude 
of correction signal is due to the attenuation by the integrator. This re- 
sult and the result from (b) show that in a closed loop analysis of the 
adaptive circuit, we need to deal with the parameter drift signal only. 

The perturbation signal does not affect the magnitude of correction signal 
even though it comes through the detector and rides on the correction signal. 
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(a) Error ,Y 2 (t) 0 

-10 
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m 



(b) Error squared 0 



time 
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(d) Correction 0 
signal YjAt) 

9 -2 




pj_g^ 12—1 Wave shapes of signals in adaptive circuit 1 with square 
wave perturbation. 
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(b) when CO^ 



IYX ENTS DIVISION OF CLEVITE corporation 




RECORDER MARK U 



20 - 



(c) when CO •> = 15 cps. 0 
* **» 



-20 




Pig. 13. Correction signal with different perturbation frequency. 
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4*3 Adaptive circuit 2 



From the simulation study of adaptive circuit it is found that a 
smoothing filter is needed to eliminate perturbation signal appearing in 
correction signal* A smoothing filter G^Cs) is inserted between and 
the integrator in adaptive circuit 1* The circuit becomes as Fig* 15(a)* 

The analog computer circuit is as Fig* 15(b) 0 This circuit is designated 
as adaptive circuit 2„ The smoothing filter will filter out completely the' 
perturbation signal appearing on the correction signal but also attenuates 
the correction signal and gives more phase lag to the correction signal* 

Too much phase lag of the correction signal tends to ms ke a closed loop 
system unstable* It causes time lag of the correcting parameter and when 
the parameter change is abrupt such as step changes excessive overshoot re- 
sults on the correction signal and sometimes it produces a sharp undershoot 
before producing correction signal in correct direction* The undershoot is 
a strong correction signal in a wrong direction* It will be studied in more 
detail later* The smoothing filter is usually chosen as a compensator for 
the adaptive loop to improve step response of the system* 

1 

(a) In the simulation* a filter of the form* — - -^7 is used* 

1 * C 2^ 



The ^2 va l^ es used are 0*2* 0*5* 1* 2* 5* 10* The same perturbation signal 
and parameter drift signal used in 4*2* are used* Taking the perturbation 
signal present in correction signal as noise* correction signal to noise ratio 
is measured at the output point of the smoothing filter G^Cs)** Fig* 15 is 
the graph of correction signal to noise ratio vs* of the smoothing filter* 
According to the graph* ^ 2 or above gives more than 10 1 signal to 

noise ratio. The perturbation waves remaining in the correction signal even 
after the filtering* will be further attenuated by the integrator following* 

°f 2 gives corner frequency of 0*5 rad* /sec* erd this gives enough 
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(a) Block diagram 




(b) Computer circuit 



Pig. 15. Adaptive Circuit 2. 
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attenuation to perturbation waves of 5 cps. present in correction signal of 
this system. Thus it is concluded that an acceptable ratio of corner fre- 
quency of smoothing filter to perturbation frequency is 1 : 60 or greater. 

But the attenuation of perturbation waves in the correction signal is not so 
important as the stability of the system, therefore, this filter is designed 
as an appropriate compensator for the adaptive loop and it naturally acts 
as a smoothing filter. 

(b) Effect of adaptive loop gain on parameter drift correction 

(1) Sinusoidal parameter drift 

The parameter is drifted sinusoidally. When the adap- 
tive loop is open the parameter drifts sinusoidally and as soon as the 
adaptive loop is connected the drifting parameter is corrected to the optimum 
value. The loop gain affects the amount of steady state error as expected. 

The necessary loop gain is decided by the required steady state accuracy. 

Fig. 17(a) shows a parameter drift when the adaptive loop is open and then 
corrected parameter condition when the loop is closed. (b) is the corres- 
ponding parameter correction signal. The loop gain is a little low and the 
steady state error is comparatively large. To reduce the steady state error, 
loop gain is increased and the parameter drift is almost completely correct- 
ed as seen in Fig. 17(c). With this high loop gain, the loop was closed 
suddenly when the parameter drift was a maximum, the parameter drift was 
corrected to the optimum value with an excessive oscillation before settling 
down to steady state; Fig. 18(a) shows the parameter drift being corrected 
and (b) is the corresponding correction signal wave shape. 

(2) Triangular Parameter Drift 

Fig. 18(c) and (d) is the triangular parameter drift being 
corrected and its correction signal. The loop gain is a little low and the 
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Pig. 17. Parameter drift correction 

(a) Parameter drift being corrected with low loop gain 

(b) Parameter correction signal 

(c) Parameter drift being corrected with a -high loop- gain 
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(a) 



(b) 




(c) 




RECORDER MARK II chart no. RA2921 32 bb 




Pig. 18. Parameter drift correction 

(a) Sinusoidal parameter drift being corrected with high gain 

(b) Correction signal 

(c) Triangular parameter drift being corrected 

(d) Correction signal 



46 



correction is not perfect. 

(3) Square wave parameter drift 

The parameter is drifted by a square wave. Fig. 19(a) 
shows the square wave parameter drift being corrected and (b) is the corres- 
ponding correcting signal, in which steady state error is small. The system 
is stable. When the loop gain is too high, the system response becomes 
oscillatory and tends to be unstable but the system will never be unstable 
because overall loop is second order system in this set-up. Fig. 19(c) 
is the parameter being corrected and (d) is the corresponding correction 
signal when the loop gain is very high. 

(c) The effect of perturbation frequency on adaptive loop step 

response . 

The most extreme case of parameter drift is step change of 
parameter due to an abrupt environmental change. Therefore, only the effect 
of perturbation frequency on adaptive loop step response is investigated. 

i 

The plant parameter is drifted by square wave of 0,01 cps. For dif- 
ferent frequencies of perturbation, maximum oveishoot and rise time are 
measured. The perturbation frequencies are varied from 1 cps. to 100 cps. 

(1) Steady state magnitude of correction signal 

In this adaptive circuit there is no attenuating com- 
ponent for the perturbation signal between perturbation signal input and the 
demodulation, therefore, the perturbation frequency does not affect the loop 
gain and the steady state magnitude of correction signal is constant with 
different perturbation frequency as seen in Fig. 13. 

(2) Maximum overshoot and undershoot 

For a very low frequency of perturbation, the correction 
signal has extremely irregular wave shapes of square waves as shown in Fig. 
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Pig. 19. Parameter drift correction 

(a) Square wave parameter drift being corrected 

(b) Corresponding correction signal 

(c) Square wa\e parameter drift being corrected with very high 
loop gain 

(d) Corresponding correction signal 
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(b) corresponding correction signal. 




20(b). The perturbation frequency is 2 cps. The correction signal has 
various amounts of overshoots and undershoots. Here, undershoot is de- 
fined as a sharp and strong correction signal to a wrong direction. The 
magnitude of overshoot and undershoot depends on the point of perturba- 
tion signal wave at which an abrupt change of parameter occurs. The ir- 
regularities of correction signal wave shapes becomes greater at lower 
perturbation frequency and with perturbation frequency below 1 cps. the 
correction signal is too irregular to have a sensible correction. The 
magnitude of overshoot and undershoot also decreases as the perturbation 
frequency increases and above 5 cps. of perturbation frequency there is 
no overshoot nor undershoot. Table 4.1 lists the measured ratio of maxi- 
mum overshoot and undershoot to steady state correction signal amplitude at 
lower perturbation frequencies. 

C0^(, cps.) Max. overshoot Max. undershoot 



1 


1.63 


4.25 


1.5 


1.28 


2.06 


2 


1.19 


1.81 


3 


no overshoot 


1.38 


4 


It 


1.18 


5 


tt 


1.03 



Table 4.1 Ratio of Maximum overshoot and undershoot to the amplitude 
of steady state correction signal. 

It is noted from Table 4.1 that the maximum undershoot ratio is much great- 
er than overshoot ratio. The sharp undershoot should be avoided, which 
makes the system unstable and leads to erratic operation. Therefore, the 
perturbation frequency should be higher than 5 cps. in this adaptive control 



system. 
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(3) Rise time 

Rise time vs. perturbation frequency is plotted as Fig. 
21. The rise time increases sharply from 0^=1 cps. to u) = 4 cps. 

Above 5 cps, rise time increases very little and stays almost constant. 

Fast rise is desirable and it is obtained at lower Co^ but it accompanies 
undesirable undershoot as seen in Table 4.1. In the following simulation 
studies, CO^ = 5 cps. is used as perturbation frequency throughout, which is 
the lowest perturbation frequency that does not accompany any overshoot or 
undershoot in correction signal. 
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4.4 Adaptive Circuit 3 

A band pass filter G^(s) was added between the squarer and the de- 
tector in adaptive circuit 2, Fig, 15. The circuit becomes as shown in 
Fig. 22, The band pass filter attenuates all unnecessary frequency compon- 
ents considerably in the signal coming from squarer and passes the components 
which contain information without attenuation. The band pass filter is es- 
sential in a multi- loop adaptive control system where different perturbation 
frequencies carry information for different parameters being controlled through 
different loops. This filter also affects the stability of the adaptive loop. 
The requirement of the band pass filter is explained in Section 3.3. To meet 
the requirement a filter of following form is used. 






G. (s) — o 2 

1 S + 2/^S +4/ 

where CO^ = perturbation frequency (center frequency) 
^ = damping ratio 

The half power frequency of this- filter is calculated as, 

H - *1 i 



(4-3) 



(4-4) 



The band width - 2 tO ^ 

To decide the jf value which gives best step response to the system, the jr 5 
is varied and the step response of the system is observed. /J 1 of 0.5 or 
0.6 gives best response. 

(a) First, only the band pass filter is simulated.^ of 0.5 is taken 
and the filter becomes as. 



U? s 



G l^ " _2 



(4-5) 



S + S + CO^ 
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G 4 (s) G 2 (s) 

Pig. 22. Adaptive circuit 3 




Pig. 23. Analog computer simulation of band pass filter G^s) 
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The analog computer simulation of the band pass filter is shown in Fig. 23. 
The frequency response and phase angle characteristic is plotted as Fig. 24. 
The' calculated half-power frequency and band width and the actual half- 
power frequency and band width are as below, 





&> H1 (rad/sec) 


4c> h2 (rad/sec) 


CO\, „ (rad/ sec) 

D • W, 


calculated 


15.7 


47.1 


31.4 


actual 


19.2 


49.5 


30.4, 



Using the result of signal analysis in Section 3.3(a), following signal 
components are expected to be present in the signal from the squarer of 
adaptive circuit 3. 



D.C 

2(0 = 1.256 rad/sec 
m 

2 to m 62.8 

CO. - CO = 30.144 rad/sec 
1 m 

(D^ + (O^ = 32.656 " 

Therefore, this filter will suppresses D.C., 2 46^, and 2 (O^ and passes 

(/).- cO , and (O .+ cO components without attenuation. 

1 m 1 m 

(b) The adaptive circuit 3 simulation was obtained by connecting 

the band pass filter simulation, Fig. 24 between M and M 2 in adaptive 

circuit 2, Fig. 15. The performance of the adaptive loop is much improved 

by the addition of the band pass filter. With proper loop gain, a square 

wave parameter drift is corrected perfectly as shown in Fig. 25, where (a) 

is the square wave parameter drift being corrected and (b) is its correction 

signal. The M is 1.08 and the rise time is 0.7 sec. Even with excessive 
pt 

loop gain, the performance is less oscillatory compared to the same case in 
adaptive circuit 2, in 4.3(b)(3). Fig. 25(c) is the parameter drift being 



55 




56 






Pig. 25. Parameter drift correction 

(a) Square wave parameter drift being corrected 

(b) Corresponding correction signal 

(c) Square wave parameter drift being corrected with excessive 

loop gain 

(d) Corresponding correction signal 
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corrected with excessive loop gain and (d) is its correction signal. 

Compare these figures with Fig, 19(c) and (d) which are the same case in 
adaptive circuit 2. Fig, 25(c) shows that the parameter drift is over- 
corrected because of excessive loop gain but the response is less oscil- 
latory. The comparison shows a definite improvement in adaptive loop per- 
formance by adding band pass filter. This filter is acting as a compensator 
in this circuit. 

(c) A comparison of demodulated outputs is made between following two 
cases in sinusoidal parameter drift. 

(1) Modulated wave is not filtered by band pass filter, (adaptive 
circuit 2). 

(2) Modulated wave is filtered by band pass filter, (adaptive 
circuit 3). 

Fig. 26(a) is for the case 1. Refer to the Fig. 22 adaptive circuit 3 
for the designation of signals at each point. For the case 1 the modulated 
signal goes into the detector without filtering by band pass filter. 

The demodulated output is Y^. The modulated signal Y^ contains many un- 
necessary signal components which do not contain any intelligence as ex- 
plained in section 3.2. Accordingly, the demodulated output contains un- 
necessary signal components. On the other hand, in the case 2, Fig. 26(b), 
the Y^ passes the band pass filter which filters out unnecessary signal 
components as explained in section 3.3. Thus, the signal Y^ contains only 
signals which contain intelligence. Note the wave shapes after filtering. 
This signal is demodulated and comes out as Y^ which has clear .cut polarity 
in accordance with the parameter drift direction. This indicates the Y^ 
for the case 2 does not contain any unnecessary signal as in the case 1. 

This verifies the signal analysis in section 3.3 and 3.4. 
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Pig. 26. Comparison of demodulated output with modulated signal filtered 
and unfiltered, 

(a) Unfiltered (b) Filtered 
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4,5 Adaptive circuit 4 

The measurement lag filter G (s) is connected between the parameter 

m 

perturbation and drift input and squarer. The detection reference phase 
lag filter G^(s) is connected between reference signal input and the de- 
tector. The G m (s) represents the phase lag in measuring error. G^(s) is 
chosen to give the same phase lag to the detection reference signal at U>^ 
as G^(s) gives to the perturbation signal (the carrier) at £ 0 ^, thus the 
carrier and the detection reference signal are in phase with each other at 
the detector The band pass filter G. (s) does not introduce any phase lag 
to the carrier at and the total phase lag introduced to carrier be- 

tween Y, and Y c , is the only phase lag due to G (s). The circuit becomes 
id. m j 

as Fig. 27 and analog computer simulation is as Fig. 28. This circuit is 

designated as adaptive circuit 4. The filters G (s) and G 0 (s) take follow- 

m J 

ing forms, 



° m <s) ■ TTTs" • 



m 



g 3 < s > -TTtTs 



(a) The process of demodulation is explained and analyzed in section 

3.4, It is verified in this simulation study. First, ^ rr 0 7 7 

*TH d J 1 . 4 

is used, which gives 45° phase lag at ^ = 31.4 rad/sec/ (5cps .) . Thus 
the carrier signal in Y c and the detection reference signal Y, are in phase 

D O 

at the detector. The modulated signal Y- is multiplied by the Y, signal 
and the demodulated output is produced. The demodulated output has one 
polarity or the other in accordance with the parameter drift direction with 
respect to the optimum parameter value. Fig. 29(a) is the modulated signal 
Yj. and (b) is the demodulated output Y^. The demodulated output has exactly 
the same wave shape as obtained by graphical multiplication of Y_ and Y, in 

j o 
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Fig. 27 Adaptive circuit 



G-mOO 
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Fig. 28 Analog computer circuit for the adaptive circuit 




(b) demodulated output. 




section 3,4 Fig. 5. This signal is integrated and develops a correction 
signal of one polarity or the other. When the polarity of detection re- 
ference signal Yg is changed, the demodulated output also changes its 
polarity. These results confirm the explanation of the demodulation in 
section 3.4. 

The modulated signal and demodulated output wave shapes for square 
wave parameter drift are shown in Fig. 30. The (a) in this figure is para- 
meter drift, (b) is the output of squarer, (c) is the modulated signal com- 
ing out of the band pass filter and (d) is the demodulated output. The 
demodulated output has one polarity or the other depending on the direction 
of parameter drift. 

(b) Next, the effect of phase difference between Y_ and Y, on the 

_> b 

demodulated output is observed. This effect is explained in section 3.4 
and the result of graphical multiplication of two signals which are out of 
phase with each other by some angle is shown in Fig. 5. The response of 
the system is observed for the following three cases, 

1. Yg and Y^ are in phase. 

2. Y fi leads Y $ by 60°. 

3. Yg lags Y 5 by 60°. 

For the case one, the parameter drift being corrected and parameter cor- 
rection signal wave shapes are as Fig. 31(a) and (b) respectively; the M 
is 1.046 and rise time is 0.9 sec. For the case 2 and 3, the results are 
the same and the result of the case 2 is shown in Fig. 31(c) and (d). The 
response is overdamped and the rise time is about 2 sec. This shows that 
the overall loop gain is reduced due to the phase difference between the 
Y c and Y, ♦ This effect is predicted in the signal analysis in section 3.4. 
The phase difference between Y^ and Y^ does not affect the system seriously 
and merely reduces the loop gain which can be easily restored. 
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MEMTS DIVISION OF CIEVITE CORfORATION CIEV 



squarer 




(c) Modulated 
signal after 
filtering 



time 



(d) Demodulated 
output 




Pig. 30. Modulation and demodulation of square wave parameter drift. 
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(a) Parameter 
drift being 
corrected 





(c) 



(d) 




RECORDER MARK H CHART NO. 




Pig. 51. Effect of phase difference between Y_ and Yg 
(a),(b) Yg and Yg are in phase. 

(c),(d) Yg leads or lags Yg by 60° . 
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4.6 Complete adaptive control circuit (adaptive circuit 5) 

Finally, a complete adaptive control circuit was simulated on an 
analog computer. The block diagram of the complete system is as Fig. 32. 
This circuit will be designated as adaptive circuit 5. The analog computer 
simulation is achieved by merely replacing the G m (s) in the adaptive circuit 
4 simulation, Fig. 28 with the model and plant simulation circuit in Fig. 7 
and the correction signal is negatively fed back to the summer, output of 
which goes to the multiplier M^. The reference signal phase lag filter 
G^(s) is not used because the carrier is expected to have phase shift of 
160° at the detector, which is nearer to 180°. The input to the plant and 
model, is R(t) = 7sinl.884t (f - 0.3 cps.). 

(a) The smoothing filter G 2 (s) is acting as a compensator. To see the 
effect of the G 2 (s) on the system response, first it was removed and the 
response of the system was observed. The adaptive loop was opened at the 
output of the integrator G^(s). The parameter correction signal developed 
with the loop open, is recorded as Fig. 33(b). Fig. 33(a) is the square 
wave parameter drift signal. To close the adaptive loop, the polarity of 
the correction signal should be reversed to have negative feedback. 

When the loop is closed, this system becomes unstable, so the loop gain 
was reduced to a very small value such that the system barely maintains a 
stable condition. The correction signal becomes as shown in Fig. 33(c). 

The magnitude of the correction signal becomes very small and the signal 
shape is drastically changed compared to the correction signal shape with 
the loop open. That is why the open loop signal analysis is not exactly 
valid when the loop is closed as mentioned in section 3. Fig. 33(d) is 
the plant output and (e) is the model output. Because of too small loop 
gain and accordingly too small correction signal, a slight correction of 
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^he complete adaptive control circuit (Adaptive circuit 









: 




Pig. 33 Signal shapes in adaptive circuit 5. 

(a) parameter drift, (b) Correction signal (open loop) 
(c) correction signal (closed loop).v (d) plant output 
(e) model output. 

l 
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parameter drift is achieved. Compare the figures (d) and (e). When the 
plant parameter changes abruptly the plant output has sharp overshoot due to 
instability of the system. 

(b) A lag filter of the form ^(s) ~ i + ^ s is inserted be- 

tween the detector and the integrator. Three different values of 2, 

5, and 10 are used and the loop gain is adjusted to give stable operation. 
With different values of the correction signal wave shapes are differ- 

ent but in all cases the critical loop gain is too low. Thus the para- 
meter drift can not be corrected completely in this circuit. Fig. 34 is 

the parameter drift being corrected and corresponding correction signal 

* 

for the three cases; Fig 34(a) is for ^ = 2, (b) is for ^ = 10, (c) is 
for '2*2 = 5. An effective compensator must be found to achieve a re- 
quired performance of the system. 
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Pig. 34 Parameter drift and correction signal 

(a)t 2 = 2 (b) r 2 = 10 



(c) 

(c) r 2 = 5 
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5. Design of a Sinusoidal Perturbation Adaptive Control System 

A design of a sinusoidal perturbation system is attempted based on the 
simulation study and using a linearization technique. The design procedure 
and considerations are virtually the same for this system as general control 
systems except some special features such as modulation and demodulation of 
error signal using perturbation signal. Finding the approximate equivalent 
transfer function of circuit components and applying linearization technique, 
a preliminary design is made on analytical basis and the preliminary design 
is simulated on analog computer and the effectiveness of the analytical 
method is checked. The final refinement of design should be done on experi- 
mental basis, 

5.1 Specification 

There are many ways to specify the performance criteria. Those speci- 
fications differ depending on the application, but at least the steady state 
requirement and dynamic or transient requirement should be specified. One 
special specification for this system is hunting loss, which is defined as 
the maximum allowable deviation of parameter caused by sinusoidal perturba- 
tion at the optimum parameter value. The choice of the hunting loss depends 
on the requirement of the system and the characteristic of the IP curve; a 
small hunting loss gives less distrubance on the output of the plant but may 
make the system lock into small local minimum in IP curve and sometimes 
does not give any information about the actual parameter drift because of 
attenuation, 

5.2 Selection of Perturbation frequency 

The perturbation frequency CU^ should be much higher than the expected 
maximum frequency of the parameter drift LU^\ this is necessary in any 
modulation system. The investigation of the effect of perturbation frequency 
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on the system response in 4.3(b) shows that a perturbation frequency of 
higher than 5 cps, does not make much difference on the system response. 

At an extremely low frequency of perturbation, the system operation tends 
to become erratic due to sharp overshoot and undershoot even though the 
rise time is very short. If the perturbation frequency is too high, the 
attenuation through the system components is severe and the higher gain of 
perturbation is required. 

5.3 Linearization of the System 

To investigate the stability and response of the system analytically, 
an equivalent signal transfer function of the adaptive loop is established. 
The representative circuit block diagram is derived in section 3.1 and 
shown in Fig. 4. The following linearization technique is applied to 
linearize the adaptive loop. In the adaptive loop, a D.C. parameter correc- 
tion signal is proportional to the plant parameter error X - X^, as shown in 
equation (3-31) section 3.4(b). This indicates a possibility of derivation 
of an equivalent linear loop. In applying the linearizing technique to any 
specific system, only the variation of the IP curve about the optimum is 
important. Therefore in the general IP curve equation, set X^ = 0 and a^ = 

0. Referring to the Fig. 27 adaptive circuit 4, which is a complete block 

t 

diagram of the system, the portion between error output and the detector is 
a carrier system. To linearize the squarer and the detector, the describing 
function or a first term of Taylor series approximation is used and to get 
an equivalent signal transfer function for band pass filter G^(s) frequency 
transformation technique is used. (Ref. 1) 

The perturbation signal does not get into the picture when a closed loop 
analysis is made as verified in simulation study 4.2. Therefore, only the 
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modulating signal (parameter drift signal) is acting in the closed loop. 

(a) The parameter drift signal acts directly on the G m (s) and no 
linearization is necessary for G m (s). 

(b) Linearization of 

Standard describing function technique is used. The input to is the 
error output Y 2 (t) and can be rewritten from equation (3-7) as, 

V 2 (t) . A p C lS in(Wje - e p ) ♦ - e m ) (5-1) 



The output from M 2 is Y^(t) equation (3-10). In this equation only the 
last term is the carrier modulated by the intelligence signal and all 
other terms does not contain any intelligence and will be suppressed by 
the band pass filter. Therefore, as far as the modulating signal is 
concerned, only the last term of the Y^(t) is the output from M 2 . The 
last term of the Y^(t) is written from equation (3-10), 



2a.A A C,X 1 sin( CO t - 0 )sin ( CO t - 0 ) (5-2) 

lmpll 1 p m m 

The equivalent multiplier gain is the ratio of the peak amplitudes of 

signal input and output terms. Thus, the peak amplitude of (5-2) is 

divided by the peak amplitudeof the modulating signal term in (5-1) and 

one obtains 



G U1 — 2a A C- 
Ml 1 P 1 



(5-3) 



This same gain also can be obtained by a Taylor series approximation. 

(c) Linearized equivalent transfer function of band pass filter G^(s), 



G^s) = 



CO* S 



s 2 + s + a ^ 
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. 



: 




The steady state phase shift due to G^(s) at S = j is 0°. 

The linearization is based on the following assumptions; 

1. Phase shift curve is linear at S = j 

2. Attenuation at S = J(A> + co ) is almost the same as at S = j u) . 

1 m 1 

Considering only the modulating signal component, the input to the band 
pass filter is the signal ( 5 - 2 ) and the output is the last term signal in 
Y^(t) equation (3-21). The output is written as, 



2a A A C X A sin( u) t - 0 )sin( CO t - 0 m - A 6 ) (5-4) 

impiij i p m m 

The input (5-2) is compared with the output (5-4), the modulating signal 

is attenuated by the factor of A^ and phase lag of A 0 is introduced by 

the Gj(s). An equivalent transfer function may now be defined such as to 

approximate the phase characteristics of the original network or, 



^1 equiv.^ 1 + ^ 



(5-5) 



1 eq. 



In this case, A is A^ which is the attenuation factor of the band pass 

filter at CO . and 'J’ is determined by choosing a representative 

1 1 • 

value of CO and 'T', is calculated such that G, . (s) gives the same 

m 1 eq. lequiv. ' ° 

phase lag at CO m as A 0 , which results in the original system. 

(d) Equivalent transfer function for M 2 . 

The detection reference signal input to M 2 is 



Y 6 (t) = BC 1 sin(^ 1 t - 0 p ) (5-6) 

and the modulated signal input to M 2 is the signal term (5-4). 

The output of M 0 is the Y_ , (t) equation (3-27) and is rewritten as, 

l /a. c. 

Y 7d.c. (t) = 1” ( 2 a i A m A p C l X l A 3 )BC l Sin( %t - 9 m - A e) ( 5 - 7 ) 

The ratio of the peak amplitude of input (5-4) to output (5-7) for M 2 is 
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(5-8) 




The components after the in the adaptive loop act upon the signal terms 
directly, so no further transforms are required. 

5.4 Stability investigation of uncompensated system and compensation 

After linearizing the required components, the block diagram of the 
system is formed with only the fixed components. All components are fixed 
except the smoothing filter G 2 (s) which is acting also as a compensator. 

The block diagram is reduced to a convenient form and the overall loop 
transfer function is obtained. The gain requirement is decided on the basis 
of static and dynamic accuracy and the gain is set at the lowest acceptable 
value. Using the preferred standard method of analysis, the uncompensated 
system is analyzed. The minimum gain permissible is that which provides 
acceptable steady state error. In most uncompensated systems this gain 
exceeds the critical gain of the system and the system becomes unstable. 
Proper compensation is necessary to make the system stable maintaining the 
required gain. Many different methods of compensation are available depend- 
ing on situations, but whenever possible, a phase lag compensation is recom- 
mended because the compensator is also acting as a smoothing filter. The 
required band width is narrow in this system because the parameter drift 
frequency is generally very low and it is effective in suppressing noise. 
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5.5 Example 

Some adaptive control systems are designed analytically and simulat- 
ed on computer. The actual performance of the system is compared with the 
one predicted by analytical method. Because of difficulty in simulation 
on an analog computer, the first example uses mild compensation and was 
simulated on analog computer. In the last example an extreme case of 
design is made and it is simulated on digital computer. 

(a) Example 1 

Design an adaptive control loop to maintain optimum for the 
following plant in which the pole is drifting with environment change. 

*V 8) = S(S + 9) 

Input to the plant; R = 4sinl.884t 

Specification; Steady state error = 10% of average velocity of the para- 
meter drift. 

Maximum overshoot = less than 20% 

Hunting loss = 10% of optimum parameter value 
Expected maximum frequency of parameter drift 
= 0.1 cps. 

Design steps 

(1) Gain requirement 

Steady state error = (Qa 

K 

v 

where; = average velocity of parameter drift 

K = static error coefficient 
v 

In this case K = 10 
v 

(2) Selection of perturbation 
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Parameter perturbation signal; 




sin UJ t 
1 



Assume 1 v of parameter adjusting signal corresponds to 10% of optimum para- 
meter value. Therefore, = 1 v 

The expected maximum frequency of parameter drift is given as Cc? m = O.lcps. 
Choose CO ^ = 5 cps. which gives tO ^ : CO ^ = 1 : 50 and gives also no sharp 
overshoot in correction signal as seen in simulation study. Thus, = 1 sin 



31. 4t 



(3) Establishing a linearized block diagram of the system 

Referring to the Fig. 4 the complete block diagram of the system, 
find each block. 

1) G m (s) is generally, from 3.1 



K 

. . . m 

G m^ S ^ = 1 +U) S 

m 



The measured phase characteristic of G^(s) is the same as the phase character- 
istic of the plant (Refer to the section 4.1(c)(2) ). Therefore, the denom- 



inator of the G (s) is the same as the plant transfer function denominator, 
m 



thus 



K 

m 

G m (s) = S(1 + 0.111S) 



I G (s) I = A = 

I m I m 



K 



m 



(j aO(i+o.nijaO 



m 



m 



A was measured as 0.086 at (j) = 0.0785 rad/sec. 

m “'m 

K is calculated as 0.122, thus 



rn 



r ( \ - 0-122 

G m'‘ S; " S(1+0.111S) 



2) M. 



The linearized Mj is obtained from equation (5-3) 



G m . = 2a.A C. 
Ml 1 p 1 
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The is assumed to be 1. The attenuation factor of the perturbation 

signal Ap was measured. The value of is different depending on the 

point of the error output signal wave at which it is measured as shown in 

Fig. 10(a). 

Maximum A = 0.0087 
P 

Minimum A = 0,00199 
P 

Average A = 0.0054 
P 

Take the maximum A^ for the first trial design, Cj = 1 v, thus 

G ffl - 0.0174 



The value of affects the overall loop gain which will be adjusted 



properly at the final design. 

3) 



G l(s) = 2 



CO s 

1 



S“ + CO^S +tt^ 



Linearized G j e q u ^ v ( s ) * s > from equation (5-5), 



G lequiv^ 1 + 



leq. ‘ 



A = A, = attenuation factor of G. (s) at CO , which is measured as 6.6 



A ^ ^G x (j lO x + j a? m ) = 2.67’ 

tan 2.67° = T', S 
leq. 

S = j CO = j0.628 
J m 

_ tan 2.67* 



‘leq. 0.628 



= 0.0603 



Thus, 



G, (s) = 
leq. ' ' 



6.6 



1 + 0.0603S 



4) M 2 

From equation (5-8), G M2 



BC, 
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B is the attenuation factor of reference signal phase lag filter G^(s). 

In this example G_(s) is not used because the G (s) gives the carrier phase 
j m 

shift of 160° which is closer to 180° phase shift. Therefore, B = 1 
= 1 v. Thus, 2 =0.5 

5) ^ 2 ^ will be a compensator in this sytem and is not decided yet. 

6) The integrator G^(s) is chosen to have following form because of diffi- 
culty in the simulation circuit as explained in section 4.2. 



K 4 

°4< s > - rns 



(4) The linearized block diagram of the system becomes as Fig. 35(a) and 
which is reduced to Figure (b) . Excluding GjCs) from the system and put- 
ting = 10, the stability of the system is investigated. To obtain K v = 
10, should be 1430. The open loop transfer function of uncompensated 
system becomes, 

G o (s > = S(1 + 10S)(1 + 0.111S)(1 + 0.0603S) (5 ' 9) 

The Bode diagram of the system is plotted as Fig. 36. The phase margin is 

zero and the system is unstable. 
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° m < s) 



M, 



G, . (s) 
lequiv 




« 4 w 



HL 



(a) 






0.007K, 



S(l+10S) (1+0.1113) (1+.0605S) 




(b) 



Pig. 35 The linearized block diagram of example 1. 
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The G^(s) also can be written in standard transient response form as, 



G o (s) = S(S + 0.1)(S + 9) (S + 16.6) (5-10) 

The root locus is also plotted as Fig. 37. 

The roots are; 

r l = 0.0336 - jO.9874 
r 2 = 0.0336 - jO.9874 
r 3 = -9.2322 

r 4 = -16.5349 (5-11) 

Two roots are located in the right half of the S plane. Thus the root 
locus plot also shows the instability of the system. 

The characteristic equation of the closed loop system is, 

S 4 + 25. 7S 3 + 151.9S 2 + 14. 9S + K = 0 (5-11) 

By Routh's criterion, the critical gain of the system is found as K c = 87.9 
The required gain of the system is 149 and exceeds the critical gain. 

A compensation is necessary to make the system stable. 

(5) Compensation 

A lag compensation is adopted, the compensator is 

1 + aTS 

G 2^ S) = 1 + TS (5- 12 ) 

Compensation 1; 

To get 5° phase margin in the final compensated system, a cross-over fre- 
quency =0.5 (phase margin is 7.5° on Bode diagram) is chosen in Fig. 

36. 

aT is calculated from the following equation, 

_1 = (5-13) 

aT 10 

83 



_1 

aT 



0.5 

10 



aT = 20 



(5-14) 



|G o (j^ c ) | = 12 db on Bode diagram Fig. 36 

This must be attenuated, therefore, 

20 log a = -12 



a = 



3.98 



Substitute this value in equation (5-14) and get T = 79.6. 

G (s) = JL+JO s 

G 2 ( } 1 + 79.6 S 



Thus, 



(5-15) 



This is the required compensator to have 5° phase margin. 

(6) Investigation of stability of the compensated system. 

Putting this G2(s) into the original system and you get the open loop 
transfer function of the compensated system as. 



r 10(1+2QS) (5-16) 

oc v } ~ S(1+79.6S)(1+10S) (1+0. 111S) (1+0. 0603S) 

The Bode diagram of the compensated system is plotted in Fig. 36. Now, 
the system has 5° phase margin and expected to be stable. 

The open loop transfer function of the compensated system is rewritten 
in standard transient response form as, 

r . v 37 . 45 (S+0. 05) (5-17) 

oc k ; " S (S+0. 0126) (S+0. 1) (S+9) (S+16.6) 

The characteristic equation of the closed loop system is 

S 5 +25.713S 4 +152.224S 3 +16.813S 2 +37.638S+1.873 = 0 (5-18) 

The roots are 

r t = -0.009526 + jO.4974 
r 2 = -0.009526 - j0.4974 
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r 3 = -0.05039 

r. = -9.0531 
4 

r 5 = -16.5905 

All roots are in the left half of the S plane. 

Appi oximated of the compensated system is found as 0.01916 

M for this = 1.9417 
Pt 

(7) Simulation of the compensated system on analog computer 

The compensated system is simulated on analog computer. The actual 
response of the system is recorded as Fig. 38 where (a) is the parameter 
drift being corrected and (b) is the corresponding correction signal. The 
response of the system is different depending on the direction of drift of 
parameter as explained in Section 4.1(b). We take the worst case of response, 
that is, the response of the system when the parameter drifts to the direc- 
tion of decreasing system . From the recorded result of the system 

response, following data are obtained. 

Steady state error = 16.57, 

M = approximately 2 

(8) Comparison of theoretical and actual data of response 



theoretical actual 



steady state 






error 


107. 


16.57. 


M ^ 


1.9417 


2.0 


Pt 
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Fig. 38 The response of the system. Example. 



(b) Example 1.1 

A different compensator is designed for the example 1 and the compensat- 
ed system is simulated on analog computer. The response of the actual system 
is also compared with the response predicted analytically. 

Compensation 2 ; 

(1) To get 10° phase margin in the final compensated system, U)^ =0.3 
(phase margin is 17° on Bode diagram) is chosen from uncompensated system 
Bode diagram. 

Using equation (5-13), aT = 33.3 

| G Q (j u> c ) | = 20.5 db (From Bode diagram of uncompensated system) 

20 log a = -20.5 a = ^ 

T = 353 

Thus, _ 1 + 33. 3S (5-19) 

G 2 C ' 1 + 353 S 

This is the required compensator to have 10° phase margin in the final 
compensated system. 

(2) Investigation of stability of the compensated system 

Putting the above G2(s) into the original system and get the open loop 
transfer function of the compensated system as, 



G (s) lo q+s ^sj 

oc* ; " S (1+353S) (1+10S) (1+0. 111S) (1+0. 0603S) 



(5-20) 



The Bode diagram of the compensated system is plotted in Fig. 39. Now, 
the system has 10° phase margin. 

The open loop transfer function is rewritten in standard transient re 
sponse form as, 



G 

oc 



(s) 



14.1 (S+0. 03) 

S (S+0. 00284) (S+0. 1) (S^ 9) (S+16. 6) 



(5-21) 
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The characteristic equation of the closed system is 

S 5 +25.70284S 4 +151.973S 3 +15.328S 2 +14. 14238+0,423 = 0 (5-22) 

The roots are 

r x = -0.02784 + j0.3025 
r 2 = -0.02784 - j0.3025 
r 3 = -0.03062 
r 4 » -9.0157 
r 5 = -16.6008 

All roots are in the left half of the S plane. 

Approximate of the compensated system is 0.0915 
M for this = 1.758 

(3) Simulation of the compensated system on analog computer 

The compensated system is simulated on analog computer. The result 
is shown in Fig. 40, where (a) is the parameter drift being corrected, 

(b) is the corresponding correction signal, (c) is the demodulated output 
and (d) is the error between the plant output and model output. 

From the results: 

Steady state error = 11.1% 




(4) Comparison of theoretical and actual data of response 



theoretical actual 



steady state 






error 


107. 


12.17. 


M * 
pt 


1.758 


2.0 



( 5 ) Simulation of compensated system on digital computer 



Because of difficulty in the simulation of the required compensator 
on analog computer, an approximate simulation of the required compensator 
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Pig. 40) Example' 1.1 response of compensated system 

(a) parameter drift being corrected 

(b) Correction signal 

(c) Demodulated output 

(d) Error .< . ■ 



was used in the simulation on analog computer in (3). 

The same system is simulated on a digital computer and the result 
is compared. Refer to the Appendix I for the FORTRAN program of the 
simulation. Fig. 41 is the graph of correction signal plotted from the 
numerical data output of the computer. If the Fig. 41 is compared with 
Fig. 40(b) there is a difference in signal shape. This is due to the ap- 
proximate simulation of the compensator in Fig. 40. 

(c) Example 1.2 

To get more phase margin in the system of Example 1, another compensa- 
tor is designed and the compensated system is simulated on a digital compu- 
ter. 

Compensation 3: 

(1) To get 40° phase margin in the final compensated system, (jq =0.1 
(phase margin is 45° on Bode diagram) is chosen from uncompensated system 
Bode diagram. Calculated in the same way as in example 1, the required com 
pensator is found as 

1 + 100S 

G 2 (s) = 1 + 10 4 S 

(2) Investigation of stability of the compensated system 

Putting this G£(s) into the original system and get the open loop 
transfer function of the compensated system as, 

G (s) = 10(1+1 . Q 0S} 

oc S(1+10 4 S) (1+10S) (1+0. Ills) (1+0.0603) 

The Bode diagram of the compensated system is plotted in Fig. 42. Now, the 
system has 40° phase margin. 

The open loop transfer function is rewritten in standard transient response 
form as 
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G (s) = 

oc -4 

S(S+10 ) (S+0. 1) (S+9) (S+16. 6) 

The characteristic equation of the closed system is 

S 5 +25.7001S 4 +151.90257S 3 +14.91519S 2 +1.49149S+0.0149 = 0 

The roots are 

r 1 = -0.04353 + jO. 08431 
r 2 = -0.04353 -jO. 08431 
r 3 = -0.01108 
r 4 = -8.9951 
r $ = -16.6069 

All roots are in the left half of the S plane. 

From the pair of complex roots, of the system is found as 0.459 and 
the M for this value of is 1.195. But this is no more true because 

these complex roots are no more dominant in the system. The dominant root 
is r^ and real. Accordingly the expected response of the system should be 
an overdamped one and the time constant is large. 

(3) The simulation of the compensated system on a digital computer 

The compensated system is simulated on a digital computer. The FORTRAN 
programming for the simulation on a digital computer is shown in Appendix I. 

In calculating for the simulation, considering an overdamped re- 
sponse, A = 0.00199 which is the minimum value of A is used and the K. 

P P 4 

is found as 6240 to get K = 10. The result of the simulation is shown in 

v 

Fig. 43 which shows the pole drift being corrected and the corresponding 
correction signal. The Fig. 44 shows the pole drift being corrected and 
the error. From the results, the steady state error can not be found ex- 
actly because the signal does not reach the steady state in 40 seconds used 
in the computer simulation. 
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{&) Comparison of theoretical and selva 1 response 

It is predicted theoretically that the compensated system response is 
©verdamped and the time constant is large® The actual response shows the 
same pattern of response as predicted® The theoretical steady state error 
is 10%, The actual steady state error could not be measured bat is i -pitted 
to reach close to the theoretical one® 

(d) Example 2 

To test the validity of the analytical method in an extreme case of 
design*, the steady state error in example 1 is changed tc 1% and rigorous 
compensation is used® All other specifications are the same as example 1® 

(1) The required K « 100 

(2) The open loop transfer function of uncompensated system is 

c m <5-»> 

o' ' " S(1 <■ 106) a •> 0.1113X1 * 0.06036) 

The Bode diagram of uncompensated system is plotted as Fig® k 5® The phase 
margin is -29° and the system is unstable® 

The G {&) is also rewritten in standard transient tespcn&e form as 9 

o 



f U3S 

^o* s ' S(S 4- 0ol){S + 5){S m 16ot) 

The characteristic equation of closed loop system is 

s 4 v 25. 7 S 3 f 151. 9S 2 *■ 14.93 ♦ 1490 = 0 

The roots are 

r 1 - 0 . 6049 4- j 2.5461 
— 0.6049 * j 2.8461 
r- * “11.209 

5 



(5*24) 



($*25) 



-15. 7008 



By Routh u s criterion, the critical gain is found as K ^ 87 0 ‘9o 

c 

The required gain 1490 exceeds the critical gain* 

(3) Compensation 

A lag compensation is adopted. To get 40° phase margin in the final 
compensated system* cross-over frequency £0 ~ 0 o l (phase margin is 45° on 
Bode diagram) is chosen from uncompensated system Bode diagram^ Fig. 45 0 

i 

Calculated in the same way as in example on® the required compensator is 
found as 8 



G 2 (s) 



1 + 100S 

1 + 10 3 s 



(5-26) 



(4) Investigation of stability of the compensated system 

Putting this (s) into the original system and get the open loop 
transfer function of the compensated system as* 



G (,) . lO Ml il O OS) (5 . 27) 

oc S(l4lo’s) (U10S> (1*0. JUS) (1+O.G6C3S) 

The Bode diagram of the compensated system is plotted in Fig. 45. how* 

the phase margin is 39 °* The open loop transfer function is rewritten in 

standard transient response form as 9 



The characteristic equation of closed system is 

sh-2 5. 70001 sh 151. 900257shl4.90151S^+1.490149S+0. 0149 - C 



(5-28) 

(5-29) 



The roots are 

r 1 = -0.04348 + j 0.08428 
r 2 « “0.04348 - JO. 08428 
= -0.01109 
r. a -8.995 

4 

“16.6069 

y} 



5 




I 



100 








All roots are in the left half of the S plane and are the same values 
as the roots in the Example 1.2. Thus the same response as in the example 
1.2 is expected, that is, the response of the compensated system is over- 
damped and the time constant is large. Furthermore, a smaller steady state 
error is expected theoretically. 

(5) Simulation of the compensated system on a digital computer 

The compensated system is simulated on a digital computer. The FORTRAN 

programming for the simulation on a digital computer is shown in Appendix 1. 

In calculating K, for the simulation, A = 0.00199 is used and the K. is 
4 P 4 

found as 62400 to get = 100. Considering the overdamped response, 
is increased by 35% of the calculated value to 84000 in order to see the 
effect of an increased gain. The system still maintains stability and the 
rise time is reduced a little compared to the response of example 1.2. The 
result of the simulation is shown in Fig. 46 which shows the pole drift 
being, corrected and the corresponding correction signal. From the result, 
the steady state error can not be found because the signal does not reach 
the steady state yet. 

( 6 ) Comparison of theoretical and actual response 

It is predicted theoretically that the compensated system response is 
overdamped and the time constant is very large. The actual response shows 
the same pattern as predicted and also shows that the 357 0 more loop gain 
than the designed loop gain does not hamper the stability of the system at 
*11, which tells that the system is in very stable condition as designed. 

The theoretical steady state error is 1%. The actual steady state error 
could not be measured but the deviation of the steady state error is ex- 
pected to be not so large. 
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6. Conclusions. 



(a) Analysis of the system . 

The open loop signal analysis showed clearly the feasibility of 
the adaptive scheme in this system. 

The signal analysis leads to the obtaining of a correct equivalent 
signal transfer block for the component in the system. 

The simulation study of the successive adaptive circuit, from the 
simplest to the complete circuit, showed the function of each component and 
the effect of component change in the circuit. 

It is confirmed that only the parameter drift signal component is 
acting through the closed adaptive loop; accordingly only the parameter 
drift signal need to be considered in the closed loop analysis. 

The simulation study showed that the analytical explanation of the 
system operations such as modulation, demodulation and filtering, are in 
close agreement with the actual operation of the system. 

(b) Design of the system . 

The obtaining of a correct error measurement transfer function 

G (s) and the attenuation factor of perturbation signal through the plant 
m 

(Ap) is most important in the analytical design because these affect the 

accuracy of the analytical design the greatest. The value of A^ differs 

depending on the point of the error output wave at which it is measured. 

When the compensated system is expected to be underdamped, high value of 

A is used in the simulation and when the compensated system is expected to 
P 

be overdamped, low value of A^ is used in the simulation for the first trial 
design. 

All the components are properly linearized and an analytical design is 
attempted. The predicted performance of the analytical design is in close 
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agreement with the actual system performance as verified in examples. 

In the case of a compensated system which has a pair of dominant 
complex roots (Example 1 and 1*1) the steady state error deviation is 
within 107 o and M deviation is within 15% between the predicted and actual 
value. 

In the case of a compensated system in which the pair of complex roots 
are not dominant and a real root is dominant, the comparison could not be 
made numerically, but the actual time performance had the same pattern as 
predicted. 

This analytical method of design simplifies the design considerably 
and is useful in obtaining a preliminary design which will be refined by 
an experimental verification. 

It is further suggested that this linearized technique is applied to 
an analysis and design of multi-loop sinusoidal perturbation adaptive con- 
trol system with some modifications for the intercoupling of adaptive loops, 
which should be further investigated. 
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APPENDIX I 



FORTRAN PROGRAMMING FOR THE SIMULATION OF 
A SINUSOIDAL PARAMETER PERTURBATION 
ADAPTIVE CONTROL SYSTEM 

To program the simulation of the system, the original block diagram 
of the system is rearranged to a proper form for programming as shown in 
Fig. 47 which is the rearranged block diagram of the complete adaptive 
control circuit. Fig. 32. Then, a simultaneous differential equations are 
set up in the form shown at the end of the Program Adapt 1 (page 109 ). 

The Program Adapt 1 in the following two pages, is the FORTRAN program for 
the simulation of the compensated system in Example 1.1. The input data 
and the graph control data are shown in proper form in page 110 . A 
sample numerical output is shown on page 111. 
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1 ifi* ^ Rearranged block diagram of a sinusoidal perturbation adaptive 
control system for the simulation on a digital computer. 
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3.250- 
3. 500 

3.750 
4.000 
4.250 
4.500 

4.750 



• 399 1 0E+O 1- 

• 39577E+0 1- 
• 30625E+01- 
. 14998E+01- 
• 38917E+00 
.21934E+01 
. 35 20 4 E+01 
•40801E+0 1 

• 375 1 2E+0 1 

• 26061 E+01 



• 381 3 1E+01 

• 39432E+01 
•32066E+01 
•17681 E+01 

• 51 772E-01 ■ 
.1 8741E+01 
•32972E+01 

• 3981 5 E+01 
•37897E+01 
•27777E+01 



. 17783E+00- 
. 14467E-01- 
•. 14408E + 00- 
■•26829E+00- 
33740E+00- 
■• 31 935E+00- 
■ .22322E+00- 
■•98508E-01- 
• 38573E-01- 
. 17165E+00- 



. 13639E+01 
. 12966E+01 
. 1 3352E+01 
. 1 35 1 5E+01 
. 14774E+01 

• 1 71 1 0E+0 1 
.17812E+01 

• 1 8987E+01 
. 18800E+01 
. 18804E+0 1 



• 381 35E+01 
•27664E+01 

• 3442 1 E+01 
•46470E+01 

• 371 3 1 E+01 
.23575E+01 
•29782E+01 
• 40989E+01 
• 33238E+01 

• 2 1 939E+01 



5.000 . 
5.250- • 

5.500- . 

5.750- . 
6 . 000 -. 
6. 250*-. 

6.500- . 

6.750- . 

7.000 . 
7.250 . 



893 16E+00 

10144E+01- 

27007E+01- 

37993E+01- 

40708E+01- 

34550E+01- 

20869E+0 1- 

26484E+00- 

16153E+01 

31434E+01 



.1 1495E+01 
• 73698E+03 
•24559E+01 
•36401 E+01 
.40252E+01 
. 35242E+01 ■ 
•22648E+01 ■ 

• 51 085E+03- 
.137 62E+01 ■ 

• 2 9605E+01 « 



• 25633E+00- 

• 27742 E+00- 

• 24480E+00- 

• 1 5925E+00- 
.45597E-01- 

• 6921 0E-01- 
*• 1 7786E+00- 

• 24601 E+00- 
■• 2391 OE+OO- 
•. 18291 E + 00- 



.201 19E+01 

• 2 1 224E+0 1 
•231 52E+01 
•23608E+01 
.23360E+01 

• 23203E+0 1 
• 22784E+01 
.240 1 3E+01 
•25276E+01 
•26517E+01 



. 27297E+01 
• 38739E+01 

• 290 1 7E +01 
.171 94E+01 
•23879E+01 
•36745E+01 

• 2951 /E+01 

• 1 685 1E+01 

• 2 1 787E+01 

• 334 1 3 E+01 



7.500 • 39867E+0 1 .38835E+01 -. 10322E+00-.27289E+01 
7.750 • 39625E+0 1 .39601 E+01 -. 2361 9E-02-.26758E+01 
8.000 • 307 54E+0 1 .31717E+01 .96302E-01-.26633E+01 

8.250 . 151 80E+0 1 .16909E+01 . 17297E+00-.26638E+01 

8.500- . 36972 E+00-. 15033E+03 .2 1 939E+00-. 27329E+01 

8.750- .21769E+01-.1 9676E+01 .20928E+OO-.28556E+O1 

9.000— .35 104E+0 1-. 33663E+01 . 14413E+00-.28980E+01 

9.250-.40788E+01-.401 72E+01 . 61 639E-0 1- .29286E+01 

9.500- . 37589E+0 1-.37887E+01 -.29821 E-01 -. 289 1 1 E+01 

9. 7 50- .262 10E+01-.27431E+01-. 1 22 07 E +00- . 28671 E+01 



•25142E+01 
. 1417 1 E+01 
• 20254E+01 
• 33272E+01 
. 25233E+01 
• 1 2438E+01 
. 177 34E+01 
• 30600E+01 
.23781 E+01 
. 1239 1 E+01 



10 . 000 - 

10.250 

10.500 

10.750 
11.000 

11.250 

1 1 .500 

11.750 
12 . 000 - 
12.250- 



• 91 222E+00- 
• 99548E+00 
• 26860E+0 1 
• 37920E+0 1 
• 40725E+01 
• 34654E+01 
• 21037E+01 
• 28433E+00 
. 1 5973E+0 1* 

• 31 309E+01- 



’• 1 0906E+01 ■ 

• 80536E+00- 

• 251 65E+01 • 

• 36806E+01 ■ 

• 404 18E+01 ■ 
•35162E+01 
•22350E+01 
.46815E+O0 

*• 14202E+01 
■•29970E+01 



■. 17834E+00- 
■•1901 1 E+00- 
16951 E+00- 
■•11141 E+00- 
-.30716E-01- 

• 50821 E-Ol- 

• 1 3 130E+00- 

• 1 8382E+00- 

• 1 771 5E+00- 

• 13390E+00- 



. 29298E+01 

• 29963E+01 
•31 074E+01 

• 31 278E+01 

• 3095 1 E+0 1 
• 30610E+01 
• 30224E+01 
•30866E+01 
•31713 E+01 
•32418E+01 



. 17242E+01 
•29897E+01 
• 2 1748E+01 
•98539E+00 
. 15418E + 01 
•29220E+01 
• 22727E+01 
. 10337E+01 
. 14486E+01 
•2738 1 E+01 
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